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(54) A windowing technique for blind equalisation 



(57) A blind convergence technique is restricted to 
using a subset of equalizer output samples. Illustrative- 
ly, a receiver implements a windowed MMA approach. 
In this windowed MMA approach, a sample window 



overlays the two-dimensional plane representing the set 
of equalizer output samples. Only those equalizer out- 
put samples appearing within the sample window are 
used during filter adaptation. 
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Description 

Background Of the Invention 

s The pres nt invention relates t communications equipment, and, more particularly, to blind equalization in a re- 

ceiver. 

In blind equalization, the equalizer of a receiver is converged without the use of a training signal. As known in the 
art, there are two techniques for blind equalization: one is referred to herein as the "reduced constellation algorithm" 
(RCA) (e.g., see Y. Sato, "A Method of Self-Recovering Equalization for Multilevel Amplitude-Modulation Systems," 

10 IEEE Trans. Commun, pp. 679-682, June 1975; and U.S. Patent No. 4,227,152, issued October 7, 1 980 to Godard); 
and the other technique is the so-called "constant modulus algorithm" (CMA) (e.g., see D. N. Godard, "Self-Recovering 
Equalization and Carrier Tracking in Two-Dimens tonal Data Communications Systems," IEEE Trans. Commun, vol. 
28, no. 11, pp. 1867-1875, Nov. 1980; and N. K. Jablon, "Joint Blind Equalization, Carrier Recovery, and Timing Re- 
covery for High-Order QAM Signal Constellations", IEEE Trans. Signal Processing, voi. 40, no, 6, pp. 1383-1398, 

is 1992.) Further, EP-A-0 806 854 presents an new blind equalization technique — the muttimodulus algorithm (MM A) 
— as an alternative to the above-mentioned RCA and CMA approaches. 

Unfortunately, whether using the RCA, CMA, or MMA, approaches, the ability to blindly converge the equalizer is 
also affected by the number of symbol levels represented in the signal point constellation. In other words, the difficulty 
of "opening the eye" (as this term is used in the art) increases when the number of symbol levels increases. 

20 

Summary of the Invention 

Any blind convergence technique is affected by the distribution of the output signals, or samples, of the equalizer. 
As such, an increase in the number of symbol levels increases the distribution of the equalizer output samples, which, 
2S in turn, makes it more difficult to blindly converge the equalizer. Therefore, and in accordance with the inventive concept, 
a blind convergence technique is restricted to using a subset of the equalizer output samples. This improves the ability 
to blindly converge the equalizer notwithstanding an increase in symbol levels. 

In an embodiment of the invention, a receiver implements a windowed MMA approach. In this windowed MMA 
approach, a sample window overlays the two-dimensional plane representing the set of equalizer output samples. Only 
30 those equalizer output samples appearing within the sample window are used during filter adaptation. Two different 
variations of the windowed MMA approach are described. 

Brief Description of the Drawing 

36 FIG. 1 is an illustrative block diagram of a portion of a communications system embodying the principles of the 

invention; 

FIG. 2 is an illustrative block diagram of a phase-splitting equalizer; 
FIG. 3 is an illustrative block diagram of a portion of an adaptive filter for use in an equalizer; 
FIG. 4 is an illustrative block diagram of a cross-coupled equalizer; 
40 FIG. 5 is an illustrative block diagram of a four-fitter equalizer; 

FIG. 6 is an illustrative signal point plot of an output signal of an equalizer before convergence; 

FIG. 7 is an illustrative signal point plot of an output signal of an equalizer for a system using the MMA blind 

equalization method; 

FIG. 8 is an illustrative signal point plot illustrating the reduced signal point constellation of the RCA blind equali- 
4S zation method; 

FIG. 9 is an illustrative signal point plot illustrating the circular contour of the CMA blind equalization method; 
FIG. 10 is an illustrative signal point plot illustrating the piecewise linear contours of the MMA blind equalization 
method; 

FIGs. 1 1 and 1 2 are illustrative block diagrams of a portion of a receiver embodying the principles of the invention; 
so FIGs. 13, 14, and 15, are illustrative signal point plots illustrating the piecewise linear contours of the MMA blind 

equalization method for a nonsquare constellation; 

FIGs. 16 and 17 are illustrative signal point plots of an output signal of an equalizer for a communications system 
using a two-step MMA blind equalization method; 

FIG. 18 shows a table providing a general comparison between the RCA, CMA, and MMA, blind equalization 
ss methods, without CHCF; 

FIG. 1 9 shows a tabl of illustrative data values for use in th RCA, CMA, and MMA, blind equalization methods; 
FIG. 20 is an illustrative graph of an incorrect diagonal solution for a 64-CAP signal point constellation; 
FIG. 21 shows an illustrative partitioning of a signal point constellation using a windowing approach; 
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FIG. 22 shows an illustrativ partitioning of a signal point constellation using a halt-constellation windowing ap- 
proach; 

FIG. 23 shows a table of illustrativ data values for use in the half-constellation WMMA approach; 
FIG. 24 shows a table of illustrative data values for comparing cost functions; 
5 FIG. 25 shows an illustrative partitioning of a signal point constellation using an dge-point constellation windowing 

approach; 

FIG. 26 is an illustrative signal point plot of an output signal of an equalizer using an LMS algorithm; and 

FIGs. 27 - 28 are illustrative signal point plots of an output signal of an equalizer using a Half-Constellation WMMA 

algorithm and Edge-Point WMMA algorithm, respectively. 

TO 

Detailed Description 

An illustrative high-level block diagram of a portion of a communications system embodying the principles of the 
invention is shown in FIG. 1 . For illustrative purposes only, it is assumed that receiver 10 receives a CAP (carriertess, 
is amplitude modulation, phase modulation) signal, which can be represented by: 

r(0=%[a n p(t-nT)-b n p(t-nT)]+W) (1) 

H 

20 

where a n and o n are discrete-valued multilevel symbols, pftland p (t) are impulse responses which form a Hilbert 
pair, T is the symbol period, and %(t) is additive noise introduced in the channel. 

It is assumed that the CAP signal in equation (1) has been distorted while propagating through communications 
channel 9 and experiences intersymbol interference (ISI). This ISI consists of intrachannel ISI (a n or b n symbols inter- 
ns fering with each other) and interchannel ISI (a„ and b n symbols interfering with each other). The purpose of receiver 
10 is to remove the ISI and minimize the effect of the additive noise \(t) to provide signal f'(f). The inventive concept 
will illustratively be described in the context of a windowed MMA blind equalization algorithm for use within receiver 
10. However, before describing the inventive concept, some background information on adaptive filters and the above- 
mentioned RCA, CMA, and MMA algorithms is presented. Also, as used herein, an adaptive filter is, e.g., a fractionally 
30 spaced linear equalizer, which is hereafter simply referred to as an FSLE equalizer or, simply, an equalizer. 

Equalizer Structures 

An illustrative phase-splitting FSLE equalizer 100 is shown in FIG. 2. It is assumed that FSLE equalizer 100 op- 
35 . erates on an input signal comprising two dimensions: an in-phase component and a quadrature component. FSLE 
equalizer 100 comprises two parallel digital adaptive filters implemented as finite impulse response (FIR) filters 110 
and 1 20. Equalizer 1 00 is called a "phase-splitting FSLE" because the two Fl R filters 1 1 0 and 1 20 converge to in-phase 
and quadrature filters. Some illustrative details of the equalizer structure are shown in FIG. 3. The two FIR filters 110 
and 1 20 share the same tapped delay line 1 1 5, which stores sequences of successive Analog-to-Digital Converter (A/ 
40 D) 125 samples r k The sampling rate f/Tbf A/D 125 is typically three to four times higher than the symbol rate 1/T 
and is chosen in such a way that it satisfies the sampling theorem for real signals.. It is assumed that T/T- i, where./ 
is an integer. 

The outputs of the two adaptive FIR filters 110 and 120 as shown in FIG. 3 are computed at the symbol rate 1/T. 
The equalizer taps and input samples can be represented by a corresponding N-dimensional vector. As such, the 
4S following relationships are now defined: 

r n T = t r *' r /r-T/ r *-wJ = vector of ND samples in delay line; (2) 

50 t 

c n = \ c o t ' c 2,> "' c n) = vector °* in-phase tap coefficients; and (3) 
d n T = [d 0 , d 1 , d 2 , d N I = vector of quadrature phase tap coefficients; (4) 

55 

where the superscript Tdenotes vector transpose, the subscript n refers to the symbol period nT, and k = in. 
Let y n and y n be the computed outputs of the in-phase and quadrature filters, respectively, and: 
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yn = C n Tf n; and < 5 ) 



5 ?« = dn T r„. . (6) 

An X/Y display of the outputs y n and y n or, equivalently, of the complex output y n + j y •„ is called a signal 
constellation. FIGs. 6 and 17 show an 64-CAP constellation before and after illustrative convergence using the MMA 
algorithm. (The term •64-CAP,' refers to the number of predefined symbols in the signal space or signal constellation 

'0 each symbol representing 6 bits since 2 s = 64. Additional information on a CAP communications system can be found 
in J. J. Werner, "Tutorial on Carrierless AM/PM - Part I - Fundamentals and Digital CAP Transmitter,' Contribution to 
ANSIX3T9.5 TP/PMD Working Group, Minneapolis, June 23, 1992.) After convergence, the signal constellation con- 
sists of a display of the complex symbols A n = a n + jb n corrupted by some small noise and ISI. 

In the normal mode of operation, decision devices (or slicers) 1 30 and 1 35 shown in FIG. 2 compare the sampled 

75 outputs y n and y m of equalizer 100 with valid symbol values a„ and b n and makes a decision on which symbols have 
been transmitted. These sliced symbols will be denoted d„ and b n The receiver then computes the following in-phase 
and quadrature errors e n and 

20 e n=Vn'^ ( 7a > 

*. (7b) 

25 and the tap coefficients of the two adaptive filters are updated using the familiar least-mean-square (LMS) algo- 

rithm, i.e., 



30 



35 



40 



45 



50 



55 



d^z-d^-a e m r m (8b) 

where a is the step size used in the tap adjustment algorithm. 

Turning now to FIG. 4, a cross-coupled FSLE, 200, is shown. For this equalizer structure, the A/D samples are 
first fed to two fixed in-phase and quadrature FIR filters, 210 and 205, respectively. In this case, the sampling rate VT 
of NO 1 25 is typically equal to four times the symbol rate 1/7. The outputs of the two fixed FIR filters are computed at 
a rate 1/7" that is consistent with the sampling theorem for analytic signals as known in the art. The output signals are 
then fed to equalizer 200 having a so-called cross-coupled structure. Typically, 1/7" is twice the symbol rate 1/7 

The cross-coupled equalizer 200 uses two adaptive FIR filters 215a and 215b, each with tap vectors c n and d,. 
For simplicity, the same tap vector notations c n and d n (which have been used for the previous described equalizer 
100 of FIG. 2) are used again. However, it should be clear to those skilled in the art that the tap vectors are different 
for the two types of equalizers. These two filters are each used twice to compute the outputs y n and y m of the equalizer. 
Let r„ and r, be the output vectors of the in-phase and quadrature filters that are used to compute the outputs of the 
cross-coupled equalizer. The following definitions can be made: 

R fl = r n +yF n ,and (9b) 

Yn-yn+fin* (9c) 

The complex output Y n of the equalizer can be written in th following compact way: 
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^ = cjR„. (io) 

where the ast risk * denot s complex conjugate. Making the following definitions for the sliced complex symbol A n and 
5 th complex error 

K=** + $* (11a) 
*»=V A n . (11b) 

The LMS algorithm for updating the complex tap vector C„ can be written as: 

C - = C -aE*FL. (12) 
n+i n n n * * 

Turning now to FIG. 5, a four-filter FSLE is shown. Four-filter equalizer 300 has the same general structure as 
cross-coupled FSLE 200 shown in FIG. 4, except that the adaptive portion consists of four different filters rather than 
two filters which are used twice. For this reason it is called a four-filter FSLE. The two output signals of equalizer 300 
are computed as follows: 

y m -clr m + d T iJt r„, and (13a) 

y n = ^lr n .dlr n . (13b) 

Using the definitions for the in-phase and quadrature errors e n and e ' n in equations (7a) and (7b), the following 
tap updating algorithms for the four filters result: 

d^-d^ + a?.^, (14b) 
02*+ 1 = c 2 „- a ?„?„, and (15a) 

di^-d^-ae,?,. (15b) 

Having generally described the structure of some prior^art equalizers as shown in FIGs. 2 - 5, a general overview 
of the concept of blind equalization will now be described using the equalizer structure of FIG. 2. 

Concept of Blind Equalization 

In the normal (steady-state) mode of operation, the decision devices in FIG. 2, i.e., slicers 130 and 135, compare 
the equalizer complex output samples, (where Y n = y n + j y n ), with all the possible transmitted complex symbols, 
A n (where A n =a n + jb n ), and selects the symbol A n which is the closest to Y n The receiver then computes an error, 
E ff where: 

ss £ n =Y«-K ( 16 > 

which is used to updat th tap coefficients of equalizer 1 00. This typ of tap adaptation is called "decision directed" , 
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because it uses the decisions of slicers 130 and 135. The most common tap updating algorithm is the LMS algorithm, 
which is a stochastic gradient algorithm that minimizes the mean square error (MSE), which is defined as: 

MSE A E[\E n \ 2 ) = E[\Y H - AJ 2 ] = E[e* n ] + E[e B >], (17) 

where E[ ] denotes expectation and e n and?, are in-phase and quadrature errors, respectively. 
At the beginning of start-up, the output signal of equalizer 100, is corrupted by a lot of intersymbol interference, 
10 as illustrated in FIG. 6. The latter represents experimental data obtained for a 64-CAP receiver using a phase-splitting 
FSLE as represented by FIG. 2. 

When a training sequence is used during start-up (i.e., a predefined sequence of A n symbols), the receiver can 
compute meaningful errors E n by using the equalizer output signal V n and the known sequence of transmitted symbols 
A n In this case, tap adaptation is said to be done with "ideal reference" to distinguish it from decision directed tap 
1$ adaptation. 

However, when no training sequence is available, equalizer 100 has to be converged blindly. In this case, a deci- 
sion-directed tap updating algorithm cannot be used to converge the equalizer, because the slicer makes too many 
wrong decisions, as should be apparent from FIG. 6. 

As such, the philosophy of blind equalization is to use a tap adaptation algorithm that minimizes a cost function 
20 that is better suited to provide initial convergence of equalizer 100 than the MSE represented by equation (17). The 
cost functions used in the RCA, CMA, and MM A algorithms are described below. 

Convergence of an equalizer during blind start-up usually consists of two main steps. First, a blind equalization 
algorithm is used to open the "eye diagram." (Hereafter, this will be referred to as "it opens the eye.") Once the eye is 
open enough, the receiver switches to a decision directed tap adaptation algorithm. 

25 

Reduced Constellation Algorithm (RCA) 

This section provides a general overview of the RCA algorithm. This general overview is then followed with a 
description of the RCA algorithm in the context of each of the illustrative equalizer structures, described above. 

30 With the RCA algorithm, the error used in the tap updating algorithm is derived with respect to a signal constellation 

that has a smaller number of points than the received constellation. As illustration, it is again assumed that the signal 
constellation comprises 64 symbols. In the RCA algorithm, the reduced constellation typically consists of four signal 
points only, as shown in FIG. 8. It should be noted that the RCA algorithm requires the use of a decision device, e.g., 
a slicer, to select the closest signal point from the reduced constellation. The error between the received sample Y n 

35 and the closest signal point A rn of the reduced constellation is the complex number: 

Er.n = * r .n = Y n - A FJt , where (18) 

40 

A r* = a rj, + jb rjl = R [sgn(^ ) + jsgn(y m )] , and (19) 

where sgn (•) is the signum function and the expression on the right corresponds to the case where the reduced con- 
stellation consists of four points. The reduced constellation algorithm minimizes the following cost function: 

45 

CF=E\\E r J 2 )-E[el+e r l,]=E[\Y n - A r J 2 ], (20) 

so where E [•] denotes expectation and where e r#n refers to the slicer error. 

Now, consider the phase-splitting equalizer structure shown in FIG. 2. Using equations (5), (6), and (20), the fol- 
lowing equations result: 

55 B r,n = Vn ~ &rs> = <Vn " *sgn(y„), (21a) 
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=)>*- K» = - /fsgn(^ ) (21b) 
The gradients of the cost function represented by equation (20) with respect to the tap vectors c n and d„ are equal to: 

V c (C^ = 2E[e,yJ l and (22a) 

V d (C/0 = 2£[e rfl rJ. (22b) 



These gradients are equal to zero when the channel is perfectly equalized, i.e. when the received samples V n are 
15 equal to the symbol values A n This condition leads to the following value of R: 

20 

For example, consider the gradient with respect to the tap vector c„ From the left of equations (21a) and (21b) 
there is the condition: £[(y n -R sgn(y n ))r n J = 0. With perfect equalization y n =a n . Also, if it is assumed that different 
symbols are uncorrelated, then: Qaj^ = k^a,, 2 ], where k„ is a fixed vector whose entries are a function of the 
channel. The above condition can then be written as: £[aj;] - R EfsgnfaJaJ = 0. Noting that sgn (a^ a n = la n' ^ 
25 solving for R, equation (23) results. 

The nonaveraged gradients in equations (22a) and (22b) can be used in a stochastic gradient algorithm to adapt 
the tap coefficients of the equalizer, so that the following tap updating algorithms result: 

30. c ^= C n- a iyn- Rs ^y f fi r n' an6 (24a). 

d^^d^aly^Rsgniy,)]^. (24b) 

55 Turning now to the cross-coupled FSLE structure illustrated by FIG. 4, the complex output Y n of this equalizer is 

computed from equation (10). Using this expression in equation (20), the gradient of the cost function with respect to 
the complex tap vector C n is: 

« V c =£[(y n ->4^rR n ]. (25) 

Assuming a perfectly equalized channel the following expression for R results: 

£i|/g 2 i E[|/g 2 i 

R - E^\\*^\-2E^\ (26) 



where the expression on the right is the same as the one in equation (23) for the usual case where E[\a n \]=E 
[\b n \]. The tap updating algorithm for the complex tap vector C„ is given by 

C (Wl =C n -a(y ( ,.^)*R n . (27) 

Turning now to the four-filter FSLE structure illustrated by FIG. 5, the outputs y n and /. of this four-filter equalizer 
structur are computed from equations (1 3a) and (1 3b). The gradients of the cost function in equation (20) with respect 
to the four tap vectors are similar to the ones given in equations (22a) and (22b) and will not be repeated here. The 
tap updating algorithms are given by: 
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(28a) 



5 



(28b) 



c 2 , n '*ly m -Rs&x{y n )]r m9 and 



(28c) 



dzn-at^-Zfsgn (>>„)]?„ 



(28d) 



where the constant R is the same as in equation (23). 

The main advantage of RCA is its low cost of implementation because it is typically the least complex of blind 
15 equalization algorithms. The tap updating algorithms represented by equations (24a), (24b), (27) and (28) are the same 
as the standard LMS algorithms represented by equations (8a) and (8b) except that the slicer uses a different number 
of points. 

The main disadvantages of RCA are its unpredictability and lack of robustness. The algorithm is known to often 
converge to so-called "wrong solutions.' These solutions are quite acceptable from a channel equalization perspective, 

20 but do not allow the receiver to recover the right data. It should be pointed out that the equalizer structure in FIG. 2 is 
much more likely to converge to wrong solutions than the structure in FIG. 4. This is due to the fact that the former has 
many more degrees of freedom than the latter. 

A wrong solution that is often observed with the equalizer structure in FIG. 2 is the so-called diagonal solution. In 
this case, the in-phase and quadrature filters both converge to the same filter, so that they both generate the same 

25 output samples. As a result, the signal constellation at the output of the equalizer consists of points clustered along a 
diagonal as illustrated in FIG. 20 for a 64-C AP signal point constellation. It has been found that frequency of occurrence 
of diagonal solutions is mostly communications channel dependent. Specifically, it is created when certain fractional 
propagation delay offsets are introduced in the channel. (As a point of contrast, FIG. 16 shows an illustrative correct 
solution for a 64-CAP signal point constellation using the MMA blind equalization algorithm.) 

30 Other wrong solutions can occur when the in-phase and quadrature filters introduce propagation delays which 

differ by an integral number of symbol periods. As an example, at a given sampling instant, a,, may appear at the output 
of the in-phase filter while appears at the output of the quadrature filter. This kind of wrong solution can generate 
points in the signal constellation at the output of the equalizer that do not correspond to transmitted symbols. For 
example, a 32-point signal constellation may be converted into a 36-point constellation and the 1 28-point constellation 

35 in FIGs. 1 3, 14, and 1 5 may be converted into a 144-point constellation. 

Constant Modulus Algorithm (CMA) 

This section provides a general overview of the CMA algorithm. This general overview is then followed with a 
40 description of the CMA algorithm in the context of each of the illustrative equalizer structures, described above. 

The CMA algorithm minimizes the dispersion of the equalized samples Y n with respect to a circle with radius R. 
This is graphically illustrated in FIG. 9. The CMA algorithm minimizes the following cost function: 



. where L is a positive integer. The case L=2 is the most commonly used in practice. The cost function in equation (29) 
is a true two-dimensional cost function which minimizes the dispersion of the equalizer complex output signal Y n with 
so respect to a circular two-dimensional contour. 

Now, consider the phase-splitting equalizer structure shown in FIG. 2. The gradients of the cost function with 
respect to the tap vectors c„ and d„ are given by: 



45 




(29) 



55 



V c (CF) = 2L X Eiq yj L - Fh)\ Y\ L ' 2 y n r n l and 



(30a) 
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10 



15 



20 
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30 



40 



V d {CF) = 2L* - R L )\r n \ L - 2 y„ rj. (30b) 

Assuming a perfectly equalized channel the following value for & results: 

where the expression on the right holds for the usual case where the statistics of the symbols a n and b n are the same. 
For L=2, the following stochastic gradient tap updating algorithms results: 

c„ + , = c„ - a {y\ + J?, 2 - R 2 )y n r H , and (32a) 

d n+ / = d„-a(^ 2 +^ 2 -^ 2 )^r n . (32b) 

Turning now to the cross-coupled FSLE structure illustrated by FIG. 4, the gradient of the cost function represented 
by equation (29) with respect to the complex tap vector C„ is equal to: 

V c (CF) = 2L X £[(| Y} L - fl^| ^ L " 2 />nl- W 
For L=2, the tap updating algorithm for the complex tap vector becomes: 

C„ +1 =C n -a(|y r J 2 l-^)/> nl (34) 



where R is given by the expression on the right in equation (31 ). 

Turning now to the four-filter FSLE structure illustrated by FIG. 5, the gradients of the cost function represented 
by equation (29) with respect to the four tap vectors are similar to the ones given by equations (30a) and (30b). For 
$ 5 L=2 % the tap updating algorithms become: 

- c/. n - <x(y 2 n + y\ - R 1 )y n r„ (35a) 

d^^d^+aiyl+yt-R 2 )?^ (35b) 

= - <*(y 2 » + Jl - r2 ) yX > and (35c) 

45 

*2.nH - d 2n - a(y 2 n + yl - R 2 ) y n r m . (35d) 

The constant Ft is the same as in equation (31 ). 

50 The main advantages of CMA are its robustness and predictability. Unlike RCA, it rarely converges to wrong so- 
lutions. For some applications, other than those considered here, it also has the advantage of being able to partially 
equalize the channel in the presence of carrier phase variations. The main disadvantage of CMA is its cost of imple- 
mentation. The CMA tap updating algorithm is more complex than that of the RCA algorithm and the MMA algorithm 
and, in addition, the CMA algorithm requires a so-called "rotator" at the output of the equalizer. As a result, once a 

55 certain degree of convergenc is achieved, the output signal of the equalizer must be count r-rotat dbefor switching 
to a decision-directed tap adaptation algorithm. Th need to use a rotator after the equalizer increas s th cost of 
implementation of CMA for some type of applications. It should be pointed out, how ver, that there are other applica- 
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tions, such as voiceband and cable modems, wher the rotator function is required anyway for other purposes, such 
as tracking frequency offset introduced in the channel. In these latter cases, the need to do a rotation does not incr ase 
the cost of impl mentation, and CMA becomes a very attractrv approach. 

5 Multimodulus Algorithm (MMA) 

The MMA algorithm minimizes the dispersion of the equalizer output samples y n and J* around piecewise linear 
in-phase and quadrature contours. For the special case of square signal constellations of the type used for 16-, 64-, 
and 256-CAP systems, the contours become straight lines. This is graphically illustrated in FIG. 10 for a 64-point 
io constellation. The multimodulus algorithm minimizes the following cost function: 

CF = E [(y l m - R^Y.)) 2 +(# - R L (Y„)n (36) 



1S 
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where L is a positive integer and R(Y n ) and ifT^take discrete positive values, which depend on the equalizer outputs 



Multimodulus Algorithm (MMA) - Square Constellations 

For square constellations, ft(YJ = i( V n )= R ~ constant, so that the cost function of equation (36) becomes: 

CF- CF, + CF Q = EUyZ - R') 1 +(# - 1?') 2 ]. (37) 



Unlike the cost function for CMA represented by equation (29), this is not a true two-dimensional cost function. 
Rather, it is the sum of two independent one-dimensional cost functions CF, and CFq. The application of the MMA 
algorithm in the context of the three illustrative types of equalizers (described above) will now be described. 

For the phase-splitting equalizer structure shown FIG. 2, the gradients of the cost function in equation (37) with 
30 respect to the tap vectors c„ and d„ are equal to: 

V c (CF) = 2L X £[( |yJ L - ft) \y\ L ' 2 y/J, and (38a) 

V d {CF) =2Lx E[Q y m \ L - lf)\y m \^y n rj. (38b) 
Assuming a perfectly equalized channel, the following value for Ft- results: 



45 The best compromise between cost and performance is achieved with i^2, in which case the tap updating algo- 

rithms become 

d m n-d m -*(yt-R z )y a r„ (40b) 



Turning now to the cross-coupled FSLE structure illustrated by FIG. 4, the gradient of the cost function represented 
& by equation (37) with resp ct to the complex tap vector C n is given by: 
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V c (CF)=2LX EI/CRJ, (41) 

where, 

5 

Ar-Kly/./^ly^J+yKI^J^/^iyj^J. (42) 

10 Assuming a perfectly equalized channel, the value for is: 

^ = _ in nJ ( (43) 
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which reduces to equation (39) for the usual case where the symbols a„and i>„have the same statistics. For b=2, the 
tap updating algorithm for the complex tap vector C n becomes: 

c n+ i= c „-«'<* R n . w 

where, 

K = (y 2 -B?)y+jXy 2 -R 2 )y- (45) 

Turning tiow to the four-filter FSLE structure illustrated by FIG. 5, the gradients of the cost function represented 
by equation (37) with respect to the four tap vectors are similar to the ones given in equations (6.5). For L=2, the tap 
updating algorithms become: 

d^-d^ + a^-* 2 )?,,^ (46b) 
c 2 .n+, = c 2 .n - a(j7j - K 2 ) y„ r„ . and (46c) 

d z „ + / = d z „ - a(yl - R 2 ) y.r. . (46d) 

The constant R is the same as in equation (39). 

45 The above-mentioned two-step blind equalization procedure utilizing the MMA algorithm is graphically illustrated 

by FIGs. 6, 7, 16, and 17 for equalizer 100. The output signal of equalizer 100, before any form of convergence, is 
shown in FIG. 6. As noted above, FIG. 6 represents experimental data obtained for a 64-CAP receiver using a phase- 
splitting FSLE as represented by FIG. 2. FIG. 7 illustrates the beginning of the MMA process convergence. As shown 
in FIG. 16, the MMA technique converges the equalizer enough to clearly illustrate the 64-symbol signal space as 64 

50 noisy clusters. Although these noisy clusters would, typically, not be acceptable for steady-state operation — the eye 
is open enough to allow the receiver to switch to a 64-point slicer and a decision-directed LMS algorithm. The end 
result is a much cleaner constellation, as shown in FIG. 17. Typically, a clean transition can be made between the two 
modes of adaptation, MMA and decision directed, when the symbol error rate is better than 10* 2 , although successful 
transitions have been observed for worse symbol error rates. It should be pointed out that the noisy clusters in FIG. 

55 1 6 could b further reduced by deer asing the step size in the MMA tap adjustment algorithm. Ind ed, in some appli- 
cations it may be possible to eliminate the switching to a decision directed tap adaptation algorithm. However, it should 
be noted that this would increase th start-up time and the required amount of digital precision. 
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The MMA algorithm for square constellations can be used without modification for nonsquare constellations. In 
this case, caution has to b xercised in the computation of the constant ft, because the discr te levels for th symbols 
a„ and b n do not all have the same probability of occurrence (described below). However, it has been found through 
computer simulations that conv rgence of the MMA algorithm is somewhat less reliable for nonsquare constellations 
s than for square constellations. This can be corrected by using the modified MMA discussed in the following section. 

Multimodulus Algorithm (MMA) - NonSquare Constellations 

The principle of the modified MMA is illustrated in FIGs. 13, 14, and 15, with respect to a 128-CAP signal constel- 
io lation. (A 128-point signal constellation is obtained in the following way. First define a 144-point signal constellation 
using the symbol levels ±1,13^,17,^,111, and then remove the four comer points in each quadrant.) Minimization 
of the dispersion of the equalizer output samples y n and y m is now done around piece wise straight lines. Again, this 
is done independently for y n and y m . The in-phase cost functions derived from equatbn (37) are: 

1$ 

CF Q = E[(yt - Rt) 2 ] if fej <K, and (47a) 

20 CF Q = El(y£ - fl£) 2 ] if |y,J > K. (47b) 

The quadrature cost functions derived from equation (37) are: 

2s CF / «£[(x^-tf) , ] »f lyj <K, and (47c) 

GF-£[(xf -**)»] if \y m \>K. (47d) 

30 The constant K\s a function of the signal constellation under consideration and is determined empirically In computer 
simulations for 1 28-CAP, a suggested value is K = 8. Two different moduli R A and ft 2 are used in equations (47) because 
the symbols a„ and b n used in the 128-point constellation have two sets of levels {±1 ,±3,±5,±7} and {±9,±11} which 
have a different probability of occurrence. More moduli can be used if there are more than two sets of symbol levels 
with different statistics. 

35 The moduli R A and R 2 in equations (47) are computed from equation (39) by evaluating the moments of the symbols 

over the set of symbol levels to which a given modulus applies (additionally described below). As an example, consider 
FIG. 1 3, which illustrates the Moduli for the in-phase dimension and which applies to the real symbols a n of a 1 28-CAP 
signal constellation. The moments of the symbols can be computed by considering the first quadrant only. Consider 
the subset of 24 symbols in this quadrant that applies to R^ . For these symbols a n = 1 , 3, 5, 7, 9, 11 ; and b n = 1 , 3, 5, 

40 7; so that each value of a n occurs with probability 4/24=1/6. Similarly, the R 2 subset has 8 symbols for which a„ =1 , 3, 
5, 7 and b n = 9, 1 1 , so that each value of a n occurs with probability 2/8=1/4. Thus, the variance of the symbols becomes: 

for ft, symbols, E[a 2 n ] = |(1 2 + 3 2 + 5 2 + 7 2 + 9 2 + 1 1 2 ) = 47.67, and (48a) 

45 

for R 2 symbols, E[a 2 ] = | (1 2 + 3 2 + 5 2 + 7 2 ) = 21 . (48b) 

Other moments of the symbols are computed in a similar fashion and then used in equation (39) to evaluate the 
50 values of the various moduli. 

The tap updating algorithms for the modified MMA algorithm are the same as the ones given in equations (40), 
(44), and (46), except that the constant R is replaced by either ft, or ft? depending on which equalizer output sample 
Y n is received. FIG. 14 illustrates the Moduli for the quadrature dimension and which applies to the symbols b n of the 
128-CAP signal constellation. It should be apparent from FIG. 15, which represents the union of FIGs. 1 3 and 14, that 
55 th in-phase and quadrature tap updating algorithms need not use th same moduli ft 1 r R% in a given symbol period. 
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Moments of Data Symbols 

The following description discusses the concept of "moments of data symbols." In particular, the closed-form ex- 
pressions for the moments E[\a n \ L ], £{lf> ft l L ]. and £[l/y L ] when the symbols a„and o„take values proportional to the 
odd integers ±1 ,±345,17, — , are presented. These expressions are then used to get closed-form expressions for the 
constants R used in the three blind equalization algorithms and illustrated in the table of FIG. 19 (described below). 

First, it is assumed that the symbols a n and b n have the same statistics, so that E[a n \ L ] = E[\b n \ L ]. Consider first 
the following known summations of powers of integers: 

£* = i/>i(m + l), (49a) 



fV = -m0» + 1)(2« + 1), (49b) 



2> 3 =4[*(m + l)] J , and (49c) 



%k* =-^iit.m + lX2m+iyCSm 2 +3m-l). (49d) 

A-l 30 

30 These summations can be used to find closed-form expressions for sums of powers of odd integers. For example, 
for power one: 



(1 + 3+5 +7) = (1 +2 + 3+4 + 5 + 6+7)-2(l + 2 + 3) 

(50) 

£(2A-l)=|>-2 £* = m 2 , 



where the two summations in the middle have been evaluated by using the closed-form expression of equation 
(49a). Similar series manipulations can be used for other sums of powers of odd integers. 

Now, consider square signal constellations which use symbols a n and b n with values ±1,±3,±5,±7,.« ±(2m-1), 
where m is the number of different symbol levels (in magnitude). As an example, for the 4-CAR 16-CAR 64-CAP, and 
4£ 256-CAP square signal constellations, m = 1, 2, 4, and 8, respectively. It is also assumed that all the symbol values 
are equiprobable. As a result, the moments of the symbols a n are: 



so 



E\\o n \) = - • X<2*-1)=«, 



(51) 



ss 



• S(2*-D 2 =^(4/n J -i), 



m 



(52) 
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E[\af]^- • £(2*-l) 3 = »*2m 2 - 1) , and (53) 



E[a 4 n ] = - • Z (2*-l) 4 =^(W-l)(12m 2 -7). (54) 
m ^ 15 

Next, consider the complex symbols A„ = a n + /£>„. Assuming that the symbols a n and b„ are uncorrected, the 
folbwing expressions for the even moments of the complex symbols result: 

Et|/»J 2 ] = 2£Ia^],and (55a) 

£I|/»j 4 ] = 2et|aj] + 2{£|a3] 2 (55b) 
20 Using equations (52) and (54) in equation (55b), results in: 

^N 4 )^* 4 ™ 2 - 1 )* 28 ™ 2 - 13 )- (56) 

25 The above results can now be used to get closed-form expressions for the constants R used in the various blind 

equalization algorithms. The following (remarkably simple) expressions for these constants result: 

EI^l 4m 2 - 1 (57) 
— 5 



35 



40 



q N 4] 26 (59) 



With respect to nonsquare signal constellations, the various symbol levels 2k - 1 for a„ and b n have a different 
probability of occurrence, even when all the complex symbols A n are equiprobable. This should be apparent from the 
1 28-point constellation illustrated by FIG. 1 5. In this case, the moments of the symbols have to be computed according 
45 to the general formula: 



E[\a/) = P l f i (.2k-l) L + P 2 f i (2k-l) L + P J f,(2k-l) L +- (60) 



where P, is the probability of occurrence of the symbol levels appearing in the corresponding summation. For typical 
32-CAP and 1 28-CAP constellations the expression in (60) is restricted to two different probabilities P 1 and 

Everything else being equal (i.e. symbol rate, shaping filters, etc.), it is possible to guarantee a constant average 
55 power at the output of a CAP transmitter if E[a n 2 - EI*?,, 2 ] = constant, independently of the type of signal constellation 
that is used. Of course, different signal const Ilations will have to use differ nt symbol values if th av rage power 
constraint has to be satisfi d. Thus, in general, a signal constellation will us symbol values X(2/r- 1 ) wh r X is chosen 
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in such a way that the average power constraint is satisfied. For simplicity, it is assumed that £{a^] = 1. F r square 
constellations, th value of X can then be determined from equation (52), to result in: 



Using this expression of k in equations (57), (58), and (59), the following expressions for the normalized constants 
10 R result: 



W jSflS (62) 



n 2 _ ,2 £[0 _ 3 12m 2 -7 
EUg 5 4m -1 



and (63) 



7b 2 _ j2 ^ M* 1 _ 156m 2 - 26 (64 . 
E[|^ 2 ] 5 4/7) 2 -1 

Similar expressions can be obtained for nonsquare constellations in a similar fashion. When the number of points 
in the signal constellation becomes very large, the following asymptotic values for the normalized constants result: 

m-»~ 1-155 7»«»-1342 7^ = 1.673. (65) 

Summary of RCA, CMA, and MMA Algorithms 

A general comparison of the RCA, CMA, and MMA techniques is shown in the table of FIG. 18. In addition, the 
35 table shown in FIG. 1 9 shows illustrative values, for signal constellations of different sizes, of the constants R, R v and 
R& which are used in the tap updating algorithms of the RCA, CMA, and MMA, blind equalization techniques described 
above. The data shown in FIG. 1 9 assumes that the symbols a n and b n take the discrete values ±1 ,±3,±5,±7,-». The 
closed-form expressions for these constants are derived as described above. 

Generally speaking, the RCA algorithm has less reliable convergence than either the CMA or MMA algorithms. 
40 As between the CMA and MMA algorithms, these algorithms have both benefits and drawbacks. For example, the 
CMA algorithm provides reliable convergence - thus avoiding incorrect diagonal solutions - but the CMA algorithm 
requires an expensive rotator. In comparison, the MMA algorithm does not require an expensive rotator but is more 
susceptible than the CMA algorithm to incorrect convergence. . 

4S Number of Symbol Levels 

Any blind convergence technique is affected by the distribution of the output signals, or samples, of the equalizer. 
As such, an increase in the number of symbol levels increases the distribution of the equalizer output samples, which 
— in turn — makes it more difficult to blindly converge the equalizer. This is illustrated by the following comparison 
so between the MMA blind equalization algorithm and the standard LMS algorithm. 

For the standard LMS algorithm, the cost function minimizes the error between the equalizer's output signals Y n 
and an unknown sequence of transmitted symbols A„ 



CF = E[(Y n -4,?)= EHy m -a m ) 2 +[(& -bj 1 ] = E[e\{LMS)^LMS)] (66) 

where 
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Y* =y*+jy„, and A m = a„ +jb H . 

In comparison, f r the MMA blind equalization algorithm, the cost function minimizes the dispersion of the con- 
5 stellation: 

CF - E Kyi -R 2 ? + (# - R 2 ) 2 ] = E[el(CF) + «i(CF)]. (67) 

where the expression of the constant R is given by: 



10 



f? = — — (68, 

Comparing the two cost functions in equations (66) and (67) it can be observed that errors have different interpre- 
tations for the LMS algorithm and the MMA algorithm. 
20 In the LMS algorithm, the error e rn (LMS) is defined as 

Grn (LMS)=y n -a n ; (69) 

25 and taps are updated in the opposite direction of the gradient: 

C n + 1 = c n " ^rJLMSJtn = c„ - u.(y„ - a„)r„. (70) 

30 When y n and a n represent the input and output of a slicer, the LMS-based error used during tap updating is equiv- 

alent to the error measured at the slicer and is a well-defined quantity. As such, an equalizer can converge to optimal 
solutions when errors are directly calculated with respect to the inputs and outputs of the slicer. 

In contrast, in the MMA algorithm the error e rn (CF) is defined differently. It should be noted that since the LMS 
algorithm uses second-order statistics of the signals, whereas MMA uses fourth-order statistics; a simplified version 

55 of the MMA algorithm, with L = 7, is used here for comparison purposes. For this one-dimensional MMA the error e rn 
(CF) of the generalized MMA becomes 



40 



45 



so 



55 



and the taps of the filter are updated as follows: 

«Wi = c n " ^rJOF)r n = c„ - n( |yj -fl)r„. (72) 

From equation (72), the taps are not exactly updated in the direction of the slicer errors. Instead, error minimization is 
done with reference to the constant ft that has statistical information about the real symbols a„. The occurrence of filter 
adaptation depends on the occurrence of fl, which depends on m. Consequently, the error e r JCF) has only a statistical 
meaning and an equalizer is not always guaranteed to converge to optimal solutions in terms of mean-square error 
(MSE). 

If a blind equalization algorithm is not optimum, then CF* 0 when an equalizer has converged. That is, residual 
values of the cost function exist. Illustratively, the MMA algorithm is examined to explore what are the residual values 
of the cost function CF 

For a blind start-up with a perfect convergence, y n -> Consequently, the cost function CF converges to CF a n ; 

CF = E[(y 2 n - F?) 2 ] -» CF Qtn = E[{a 2 n - F?)\ (73) 
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This cost function CF &Ji is expanded and simplift d as follows: 

CF.^EKal-R 2 )*] (74a) 

= £K 4 -2a^ 2 + * 4 ] (74b) 

-*&:-2«£|H*}+* 4 ] (74c) 

SKI 

= E[aU-2E[al)$£{+R' (74d) 

= A 4 -£[o 4 ] (74e) 

= |K|-£[a:]. (740 



It should be noted that only analysis (or the in-phase dimension is provided, and that the same analysis applies 
20 to the quadrature phase dimensbn. The cost function can be expressed as a function of the number of symbol levels 
m. From the above description of the calculations of "Moments of Data Symbols,' the constant R is given by: 



25 



fl^^ 2 . (75) 



and the fourth moment of the symbols a„ is given by: 



30 1 " J 15 



Then the cost function CF a n can be rewritten as: 



£Ia^] = -^(4/n 2 -1)(12m 2 -7). (76) 
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40 



= ( 12/W S " 7 ) 2 -^(4m 2 -l)(12m 2 -7) (78) 

= %2/n 2 -7)(m 2 -l) (79) 
75 

= 16(12/n 4 -lW + 7). (80) 
75 

<5 Equation (80) gives a simple way to express the cost function after convergence and steady-state values of the 

cost function CF a n can be easily calculated. The number of symbol levels m (in magnitude) can be computed from the 
number of constellation points Cfor C-CAP: 

so m = ^p. (81) 

Equation (80) shows that CF^ n = Ofor m=1 t and CF an * Ofor m * 1. For instance, calculating CF^ n gives the 
following results: CF^ n = Ofor 4-CAP with m= 7, CF^ n = 14.2 dB for 16-CAP with m=2, and CF^ n = 27.7 dB for 64-CAP 
55 with m=4, etc. This means that the optimum conv rgence for a blind equalizer can only be achieved for 4-CAP with 
m= 1. Residual values of the cost function CF^ n significantly increas with increasing m. Ultimately, due to large values 
of the number m, residual values of CF a ^ become so large that a blind equalizer fails to converg . 

Residual values of CF a n are increasing functions of the number m, and convergence of an equalizer is directry 
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affected by those values. In conclusion, the reliability of a blind algorithm is highly d graded with incr asing values f 
m. When th residual values f CF a n increas beyond some quantity, th eye diagram fails to open. It has been 
experimentally found that a standard MMA is only effective for CAP applications with m less that eight which corre- 
sponds to 256-CAP. 

5 

A Windowing Approach to Blind Convergence 

In accordance with the inventive concept, a blind convergence technique is restricted to using a subset of the 
equalizer output samples. This improves the ability to blindly converge the equalizer notwithstanding an increase in 
10 symbol levels. 

In an embodiment of the invention, a receiver implements a windowed MMA (WMMA) approach. In this WMMA 
approach, a sample window overlays the two-dimensional plane representing the set of equalizer output samples. Only 
those equalizer output samples appearing within the sample window are used during filter adaptation. An example is 
shown in FIG. 21, where a sample window is defined by two dotted lines along each dimension — dotted lines 601 

is and 602 for the in-phase dimension, and dotted lines 603 and 604 for the quadrature phase dimension. These dotted 
lines form both an exclusion area 600 and a sample window within the signal space. With WMMA, only those samples, 
y n w , falling within the window, i.e., outside the exclusion area, are used during filter adaptation. This is in contrast to 
MMA, which uses all of the samples. Thus, the size of the window determines the set of data, y nw , used during filter 
adaptation, and therefore effects the convergence of the equalizer. 

20 As illustration, two different variations of the windowed MMA approach are described below. The first variation is 

"Half-Constellation WMMA" and the second variation is "Edge-Point WMMA." For the purposes of these examples, 
the size of the window is varied as a function of the value of the single parameter m w which results in a square exclusion 
area. However, it should be noted that each dotted line could be associated with a different parameter value thus 
yielding non-square exclusion areas. 

25 

Half-Constellation WMMA 

A half-constellation window is shown in FIG. 22 for a 64-CAP constellation. In the case of the in-phase dimension 
the samples y„are divided into two sets by the window boundaries m w with ly n l £ rn^and ly„J > m w With new samples 
30 y n w the cost function CF is redefined as: 

CF=E[(y 2 nw -f? w f]. (82) 

35 Note that the constant R is changed to R w because the signals y nw converge to a different constellation with the 

symbols a nw where a n w = {±5, ±7}. The taps of he equalizer are updated with samples y n w and remain unchanged 
with samples ly n l. For the half-constellation WMMA, the window boundary m w is defined as: 

40 m w =m t (83) 

where m indicates the number of symbol levels, and the magnitude of the highest symbol level is given by 2m- 
1. The window boundary m w is defined in such a way that the same number of inner-point and outer-point symbols 
are included around the constant R. In other words, the data, y nw , used to update the taps are symmetrically distributed 
45 on both sides of R. It is called half -constellation WMMA because the number of partial symbols, a nwi that construct 
the new constellation is half the number of the original symbols For the half -constellation WMMA, the constant R w 
needs to be evaluated with respect to symbols a nw . With the samples, y n , w the cost function CF W now converges to 
symbols a n v ^ 

CF = E[(y 2 w - F? w ) 2 ] -> CF = E[(a 2 nw - fl^ft. (84) 
Then the constant R w is computed as 

55 
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«*=^4^ (85) 



Note that the initial index of a n w does not start with one. Therefore, the moments for the symbols a n w must be 
capable of derivation with an arbitrary initial index. The following example is for the calculation of £[a„ 2 J- The second- 
order expectation is rewritten as: 

10 

15 The parameter w denotes the number of symbol levels required for the constellation including a n vr For the half- 

constellation WMMA, 



w = lm. (87) 

Equation (87) means that half the number of the original symbol levels m is required. The parameter M denotes 
the number of the largest symbol levels and the parameter M w denotes the numberof symbol levels below the window 
boundary m w The parameters are given by: 

M=2m • 1 and M w - m - 1. (88) 

The constant R w is then calculated as: 



ft 



(89b) 
(89c) 
(89d) 



ft (48m 4 - 40/h 2 + 7)-$ (3m* -10m 2 + 7) 
?(4m 2 -l)-?(m 2 -l) 

45 

?(7m 2 -l) 

93m' -70m 2 +7 
35m 2 -5 

Values of the constant for MMA and the constant fl w for WMMA are provided in the Table shown in FIG. 23. It 
can be observed from the Table of FIG. 23 that the values of R w are always larger than those of R. 
55 Next, the residual values of the cost function for the half-constellation WMMA are computed to show how much 
reduction can b obtained. From equation (74e): 
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CF a/ , = ft 4 - - 0F m = R 4 W - E[a*J. (90) 

Replacing fl£ with: 



10 and replacing E^e^ J with: 



^ i = 93m 4 -W+7 (91) 
35/n 2 -5 



gl<J= 2(93, " 4 ^ 0, " 2+7) . < 92 > 



the cost function CF an ioT the half-constellation WMMA is then obtained as 



= -£[<„] (93b) 
_ (93m' -70m 2 + 7) 2 _ 2 # 93m 4 -70m 2 +7 (93(;) 



(35m 2 -5) 2 m 

_2_ # (93m 4 - 70m 2 + 7)(m + 2)(m - 2X1 7m 2 - 2) 
75 * 7m 2 -1 



(93d) 



A comparison of residual values of cost functions for MMA and half-constellation WMMA is provided in the Table 
shown in FIG. 24. This table shows that the cost function for 16-CAP becomes zero, and that a reduction of about 5 
dB can be obtained for other CAP systems. Thus, by using WMMA, the cost function CF gn becomes optimum for 
1 6-CAP. For other CAP systems, the cost function CF ajn is reduced. Reduction of the residual values of the cost function 
3$ CF g/ j leads to improved reliability and convergence rate of blind equalizers. 

Edoe-Point WMMA 

Edge-point WMMA is proposed as the second application of WMMA. The cost function for half -constellation WMMA 
40 in equation (81 ) can be directly applied to edge-point WMMA, except that a modification is made for the sample window 
parameters. The window parameters are illustrated in FIG. 25, where the window boundary rn^is defined as: 



m„ = 2(m-1). (94) 

With such a definition, the symbols a nw are given by a 2m . v The symbols used for filter adaptation are only those 
that have the largest values. FIG. 25 shows that these symbols are geometrically located at the edge of the original 
constellation. Because only one symbol level is involved, the calculation of R w is simply given by: 

The above equation yields the following equality: 

"fl£=e[a*J. (96) 
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Equation (96) further leads to th result: 



^.^eKa^-^l^O. (97) 

5 

Equation (97) shows that zero value can be achieved with such a cost function. That is, with edge-point WMMA, 
the cost function becomes optimum for any CAP system. 

The edge-point and half -constellation MMA are basically the same except that they use different window param- 
eters. However, the difference in parameters results in different performance. Theoretically, optimum convergence can 
10 only be achieved for 1 6-CAP with half-constellation WMMA, and can be achieved for any CAP application with edge- 
point WMMA under conditions to be discussed later 

The design parameters of the edge-point constellation are simple and easy to implement. However, expected 
performance cannot be achieved for high level CAP applications because other factors also affect convergence, such 
as the lack of enough data samples y n 

15 

Filter Adaptation 

This section describes the algorithms for updating the tap coefficients for the Half-Constellation WMMA and the 
Edge-Point WMMA approaches. For simplicity, the previous analysis for WMMA was only given for the in-phase di- 
20 mension. The complete two-dimensional cost function for WMMA is given by: 



CF = [(yl w -Rl) 2 < 98 > 

25 The cost functions can be applied to both half-constellation and edge-point WMMAs by using different definitions 

for y n ^ as described above. The gradients of the cost function in equation (98) with respect to the tap vectors c„ and 
d n are equal to: 



30 



v.- Ok- v i-C52--<)5Ui- (99) 

The taps of the filter are then updated in a stochastic fashion in the opposite direction of the gradient: 

35 f&y«j* (100) 



d^-d.-nGS.-O?^.. < l01 > 

Note that y n w and R w are different for the two versions of WMMA. 

In the implementation of the algorithms, the samples y aMf are normally calculated by using a comparator, or a look- 
up table. Alternatively, a nonlinear function /Q is proposed to determine partial samples y nvr The mction f() is defined as: 

ffYn) = J[1 + sgn(f n - m 2 J), and (1 02) 

- + sgrAyl -«;)]■ (i 03) 



40 



45 



So that: 
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/(>-,) = 



0if W^ Ay)J° if W 5 "- (1 04) 

1 otherwise • /v - r " ; I I otherwise * • 



where m w = mfor half -constellation WMMA, and m w = 2(m-1) tor edge-point WMMA. The use of the nonlinear 
equation ((■) gives the following relation: . 

,o /(>'-) = >'..-; and f(y„) = y H . w . (105) 

The cost function CF can be rewritten as: 

CF = C/O'-Xyi - Bl) 2 + /(£.)(# - £) 2 J- (106) 

and the corresponding tap updating algorithm becomes: 

°«1 = <=„ - W>S - ^)V„. (107) 



<*~, = d. - n/(£,)(?; - £)>yr. . (108) 



In the case of the in-phase dimension, both equations (100) and (107) can be used to do filter adaptation. 
FIGs. 16, and 26 - 28, show illustrative plots of converged signal constellations for the various algorithms based 
on computer simulation. FIG. 16 shows the signal constellation after converging with MMA, and FIG. 26 shows the 

30 signal constellation after converging with LMS. These two figures together show that even though initial convergence 
can be achieved with blind algorithms, the LMS algorithm is needed to obtain optimal convergence. FIG. 27 shows the 
converged constellation with the half-constellation WMMA, whereas FIG. 28 shows the converged constellation with 
the edge-point WMMA. From FIGs. 27 and 28, it can be observed that the convergence performance is improved by 
using the half-constellation WMMA, and is further improved by using the edge-point WMMA. In fact, a comparison of 

35 FIGs. 26 and 28 shows that almost the same convergence performance is achieved by using the edge-point WMMA 
and the LMS algorithm. It should be noted that in this example the step size used in the edge-point WMMA is about 
five times larger than for the other algorithms. A step size of u. =0. 00001 is used for MMA. 

It is recommended that the use of windowed MMA be limited to applications with a limited number of symbol levels. 
For very large constellation, a good performance may be difficult to obtain because of the lack of sufficient equalizer 

40 output samples during filter adaptation. 

Illustrative embodiments of the inventive concept are shown in FIGs. 11 and 1 2. FIG. 1 1 illustrates an embodiment 
representative of a digital signal processor 400 that is programmed to implement an FSLE in accordance with the 
principles of the invention. Digital signal processor 400 comprises a central processing unit (processor) 405 and mem- 
ory 410. A portion of memory 410 is used to store program instructions that, when executed by processor 405, imple- 

4 $ ment the windowed-MMA type algorithm. This portion of memory is shown as 411 . Another portion of memory, 412, is 
used to store tap coefficient values that are updated by processor 405 in accordance with the inventive concept. It is 
assumed that a received signal 404 is applied to processor 405, which equalizes this signal in accordance with the 
inventive concept to provide a output signal 406. For the purposes of example only, it is assumed that output signal 
406 represents a sequence of output samples of an equalizer. (As known in the art, a digital signal processor may, 

£0 additionally, further process received signal 404 before deriving output signal 406.) An illustrative software program is 
not described herein since, after learning of the windowed-MMA type algorithms as described herein, such a program 
is within the capability of one skilled in the art. Also, it should be noted that any equalizer structures, such as those 
described earlier, can be implemented by digital signal processor 400 in accordance with the inventive concept. 
FIG. 12 illustrates another alternative embodiment of the inventive concept. Circuitry 500 comprises a central 

55 processing unit (processor) 505, and an equalizer 510. The latter is illustratively assumed to b a phase-splitting FSLE 
as described above. It is assumed that equaliz r 510 includes at least one tap-coefficient register for storing values 
for corresponding tap coefficient vectors ( .g. , as shown in Fl G. 3). Processor 505 includes memory, not shown, similar 
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to memory 410 f FIG. 11 for implementing the windowed-MMA type algorithms. Equalizer output signal 511 , which 
represents as quenceof qualiz r output samples, is applied to processor 505. Th latt ranalyz s equalizer input 
signal 51 1 , in accordance with the inventive concept, to adapt values t the tap coefficients in such a way as t converge 
to a correct solution. 

5 The foregoing merely illustrates the principles of th invention and it will thus be appreciated that those skilled in 

the art will be able to devise numerous alternative arrangements which are within its scope. 

For example, although the invention is illustrated herein as being implemented with discrete functional building 
blocks, e.g., an equalizer, etc., the functions of any one or more of those building blocks can be carried out using one 
or more appropriate programmed processors. 

10 in addition, although the inventive concept was described in the context of an FSLE, the inventive concept is 

applicable to other forms of adaptive filters, such as, but not limited to, a decision feedback equalizer (DFE). The 
inventive concept is applicable to all forms of communications systems, e.g., broadcast networks, e.g., high-definition 
television (HDTV), point-to-multipoint Networks like fiber to the curb (mentioned above), signal identification, or clas- 
sification, applications like wire-tapping, etc. 

is Also, although the inventive concept was described in the context of a modified MM A algorithm, the inventive 
concept is applicable to other forms of equalization. 
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1. A improved method for performing equalization in a receiver, the improvement comprising: 
using a subset of output samples from an equalizer for converging the equalizer. 

2. The method of claim 1 wherein the using step includes the steps of: 

partitioning a signal point space into an exclusion region and a sample window region; and 

using that subset of output samples from the equalizer that fall within the sample window for converging the 

equalizer. 

30 3. The method of claim 2 wherein the partitioning step includes the step of 

selecting the exclusion region and the sample window region such that along one dimension of the signal 
point space a number of symbols within the sample window region is equal to one half of a number of symbols in 
that dimension. 

3S 4. The method of claim 2 wherein the partitioning step includes the step of 

selecting the exclusion region and the sample window region such that the sample window region only in- 
cludes outer-most symbols of the signal point space. 

5. The method of claim 1 wherein the using step includes the step of: 

40 adapting a set of tap coefficients of the equalizer in accordance with a multimodulus-based blind equalization 

technique operative on the subset of output samples. 

6. The method of claim 1 wherein the using step includes the step of: 

adapting a set of tap coefficients of the equalizer in accordance with a reduced constellation-based blind 
45 equalization technique operative on the subset of output samples. 

7. The method of claim 1 wherein the using step includes the step of: 

adapting a set of tap coefficients of the equalizer in accordance with a constant modulus-based blind equal- 
ization technique operative on the subset of output samples. 



so 



8. Apparatus for use in performing blind equalization in a receiver, the apparatus comprising: 



a memory for storing a convergence algorithm and for storing a set of tap coefficient values; and 
a processor a) for filtering an input signal as a function of the set of stored tap coefficient values to provide a 
55 plurality of output samples, and b) for executing the convergenc algorithm to adapt th set of stor d tap 

coefficient values as a function of a subset of th plurality of output samples. 

9. The apparatus of claim 8 wherein the plurality of output samples hav associated coordinate values within a signal 
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point spac and the subset of output samples are those utput samples having coordinate values outsid f an 
exclusion region of the signal point space. 

10. The apparatus of claim 8 wherein th convergence algorithm is a muttimodulus-based algorithm. 

5 

11. The apparatus of claim 8 wherein the convergence algorithm is a reduced constellation-based algorithm. 

12. The apparatus of claim 8 wherein the convergence algorithm is a constant modulus-based algorithm. 

10 1 3. The apparatus of claim 9 wherein the exclusion region is such that along one dimension of the signal point space 
a number of symbols within the exclusion region is equal to one half of a number of symbols in that dimension. 

14. The apparatus of claim 9 wherein the exclusion region does not include outer-most symbols of the signal point 
space. 

15 

15. The apparatus of claim 8 wherein the processor adapts the set of tap coefficients as a function of a form of blind 
equalization algorithm. 

20 
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FIG. 3 
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AFTER CONVERGENCE WITH HALF-CONSTELLATION WMWA 



6- 



-10 





• 

# 4 * if 




> 


n 


% V II 








* * * tf 


J* * 




* 


<fe * 


* * 


i 












* * * 


* #• 






* <* *. 














0 


I 1 
-5 0 


1 

5 


i 



REAL 



IMAGINARY 



10 



8- 
6- 
4 

2- 
0- 

-2- 
-4- 

-6 - 
-8 - 
-10 



FIG. 28 

AFTER CONVERGENCE WITH EDGE-POINT WMMA 



j* • n * * * * » •* 

* * .* *• 

* 4 * * *. ,ik * * 

* * X * « * * : * 

* i? * * * * £ .* 

•It * » * *f » ff 



-5 



— I — 

0 
REAL 



"T" 

5 



10 



38 



Europalsches Patentamt 
European Patent Offlc 
Office europeen des brevets 



iniiiiwiiiiiiiii 

(id EP 0 845 891 A3 



(12) EUROPEAN PATENT APPLICATION 

(88) Date of publication A3: (51) Intel. 7 : H04L 25/03, H04L 25/49 

21.03.2001 Bulletin 2001/12 

(43) Date of publication A2: 

03.06.1998 Bulletin 1998/23 

(21) Application number: 97309276.0 

(22) Date of filing: 18.11.1997 



(84) Designated Contracting States: 

AT BE CH DE DK ES Fl FR GB GR IE IT LI LU MC 
NLPTSE 

Designated Extension States: 
AL LT LV MK RO SI 

(30) Priority: 27.11.1996 US 757207 

(71) Applicant: LUCENT TECHNOLOGIES INC. 
Murray Hill, New Jersey 07974-0636 (US) 



(54) 

(57) A blind convergence technique is restricted to 
using a subset of equalizer output samples. Illustrative- 
ly, a receiver implements a windowed MMA approach. 
In this windowed MMA approach, a sample window 



(72) Inventors: 

• Werner, Jean-Jacques 
Holmdel, New Jersey 07733 (JP) 

• Yang, Jian 

Marlboro, New Jersey 07746 (US) 

(74) Representative: 

Watts, Christopher Malcolm Kelway, Dr. et al 

Lucent Technologies (UK) Ltd, 

5 Mornington Road 

Woodford Green Essex, IG8 0TU (GB) 



overlays the two-dimensional plane representing the set 
of equalizer output samples. Only those equalizer out- 
put samples appearing within the sample window are 
used during filter adaptation. 



A windowing technique for blind equalisation 



FIG. 21 



CO 

< 

00 
CO 

o 

Q. 
HI 




Pitted by Jouw. 76001 PARIS (FR) 



EP 0 845 891 A3 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 97 36 9276 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to ctaim 



CLASSIFICATION OF THE 
APPLICATION (IntC! .6) 



line 20 
line 20 
line 46 



line 27 * 
line 33 * 
column 2, 



line 57 - column 2, 
line 20 - line 25 * 



JABLON N K: "JOINT BLIND EQUALIZATION, 
CARRIER RECOVERY, AND TIMING RECOVERY FOR 
HIGH-ORDER QAM SIGNAL CONSTELLATIONS" 
IEEE TRANSACTIONS ON SIGNAL 
PROCESSING, US»IEEE, INC* NEW YORK, 
vol. 40, no. 6, 1 June 1992 (1992-G6-G1), 
pages 1383-1397, XP080365227 
ISSN: 1053-587X 

* abstract * 

* page 1383, column 2, 

* page 1384, col inn 1, 

* page 1384, column 1, 
line 32; figure 4 * 

* page 1385, col win 1, 
line 55 * 

* page 1388, column 1, 

* paragraph [CONCLUSIONS] 

SETHARES W ET AL: "Approaches to blind 
equalization of signals with multiple 
modulus" 

ICASSP-89: 1989 INTERNATIONAL CONFERENCE 
ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING 
(IEEE CAT. N0.89CH2673-2), GLASGOW, UK, 
23-26 MAY 1989, 

pages 972-975 vol.2, XP002154011 
1989, New York, NY, USA, IEEE, USA 

* abstract * 

* page 972, 

* page 972, 
line 8 * 

* page 973, 
page 974, 



column 1, line 22 
column 1, line 43 



line 30 * 
column 2, 



column 2, line 14 - line 44 * 

column 1, line 1 - line 12 * 

* page 974, column 2, line 50 - page 975, 
column 1, line 19 * 



The present search report has been drawn up for all claims 



I- 4,6-9, 

II- 15 



H84L25/03 
H04L25/49 



1.5,7,8, 
18,12,15 



TECHNICAL FIELDS 
SEARCHED (lnt.CL6) 



H04L 



| 

at 

8 



o 
8 



BERLIN 



Dale of completion of ihe seerdt 

19 January 2081 



Binger, B 



CATEGORY OF CITED DOCUMENTS 

X : particutany relevant it taken alone 

Y : particularly relevant a combined with another 

document of the same category 
A : teennotoflpea) background 
O : nofMtfrittBn (fis closure 
P : intermediate (tocumera 



T : theory or principle underlying the invention 
E : earlier patent document but pubGshed or, or 

after the filing date 
D ; document cted in the application 
L : document cited for other reasons 



& : member ol the same patent family, corresponding 
document 



2 



EP0 845 891 A3 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 97 30 9276 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION {lntCt.6} 



AXFORD R A JR ET AL: "The transversal and 
lattice multiple modulus algorithms for 
blind equalization of QAM signals" 
MILCOM 95. UNIVERSAL COMMUNICATIONS. 
CONFERENCE RECORD (CAT. NO.95CH35750) , 
PROCEEDINGS OF MILCOM x 95. SAN DIEGO, CA, 
USA, 5-8 NOV. 1995, 

2, XP902154812 
USA, IEEE, USA 



1,5,7,8, 
10,12,15 



pages 585-591 vol 
1995. New York, NY, 
ISBN: 9-7803-2489-7 

* abstract * 

* paragraph [SETUP] * 

* page 587, column 1, line 21 
column 2, 
column 1, 
column 1, 



* page 587, 

* page 588, 



page 589, 
line 7 * 



line 16 
line 34 
line 43 



line 30 * 
line 23 * 
line 53 * 
column 2, 



CHEN S ET AL: "MULTI-STAGE BLIND 

CLUSTERING EQUALISER 0 

IEEE TRANSACTIONS ON 

COMMUNICATIONS, US, IEEE INC. NEW YORK, 

vol. 43, no. 2/G4, PART 02, 

1 February 1995 (1995-02-01) , pages 

701-705, XPO00502573 

ISSN: 0090-6778 

* abstract * 

* page 702, column 2, line 50 - page 703, 
column 2, line 4 * 



1.2.7-9J 
12,15 



TECHNICAL FIELDS 
SEARCHED (IntCU) 



The present search report has been drawn up for all claims 



BERLIN 



Date d coaptation erf Uio soardi 

19 January 2001 



Binger, B 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant I taken atone 

Y : particularly relevant I combined with anottwr 

document d the same category 
A : technological background 
O : non-written disclosure 
P : tntBimedialB document 



T : theory or principle undertyftg (he invention 
E : eariar patent document, but published on, or 

after the fifing date 
D : document cited in the applcation 
L : document cited for other reasons 

& : rnernber of the same patent famOy, correspondJnQ 
document 



3 



